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SECTION  I 


INTRODUCTION 

The  problem  of  general  instability  of  a flat,  rectangular  sandwich 
panel  subjected  to  combined  inplane  edge  loads  has  been  considered  in 
Reference  1.  This  report  describes  a computer  program  which  implements 
that  analysis.  The  program  is  suitable  for  interactive  or  batch  processing, 
and  is  readily  adapted  for  use  in  interactive  design/optimization  programs. 

1.  1 PROBLEM  DESCRIPTION 

The  sandwich  panel  under  consideration  is  flat,  rectangular  and 
composed  of  three  layers  of  uniform  thickness  (Figure  1).  Each  edge  of 
the  panel  is  simply  supported,  as  indicated  in  Figure  2,  and  the  transverse 
shear  strains  of  the  core  parallel  to  the  boundary  are  assumed  to  vanish  on 
each  edge.  Edgewise  loads  applied  to  the  panel  consist  of  compression, 
shear  and  bending  forces  (Figure  3). 

Each  face  sheet  of  the  sandwich  is  idealized  as  an  isotropic  thin 
plate  which  deforms  according  to  the  Love-Kirchhoff  assumptions.  The 
thickness  and  material  properties  of  the  two  face  sheets  are  independent  of 
one  another. 

The  sandwich  core  is  assumed  to  be  incompressible  in  the  direction 
normal  to  the  plane  of  the  panel,  and  rigidly  bonded  to  the  face  sheets. 

Since  the  inplane  extensional  and  flexural  stiffnesses  of  most  core  materials 
are  generally  small  in  comparison  with  the  extensional  stiffness  of  the  face 
sheets,  the  core  strain  energy  is  taken  to  consist  solely  of  contributions 
due  to  transverse  shear  deformation.  The  elastic  properties  of  the  core 
therefore  consist  of  the  two  independent  transverse  shear  rigidities. 
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Figure  2.  Simply  Supported  Boundary  Condition. 


1.  2 METHOD  OF  SOLUTION 


The  analysis  contained  in  the  subject  computer  program  is  based 
upon  a discrete  form  of  the  principle  of  minimum  potential  energy.  By 
relating  the  components  of  displacement  in  the  face  sheets  to  the  displace- 
ments within  the  core,  the  potential  energy  associated  with  the  buckling 
deformation  of  the  panel  is  expressed  as: 

it  = it  |0  (x,y),  ij)  (x,y),  w(x,  y ),  e , e "I  (1.1) 

ppl  x y 


Here  0 (x,  y)  and  ifr  (x,  y)  are  related  to  rotations  within  the  core,  and  w(x,y) 
is  the  transverse  displacement.  Two  additional  parameters  e^  and  e^ 
describe  the  effects  of  unbalanced  face  sheets  upon  the  buckling  displace- 
ment pattern.  Equation  1.  1 is  cast  into  a discrete  form  by  the  introduction 
of  assumed  displacement  modes  satisfying  the  conditions  of  simple  support: 


0 (x,  y) 
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£ £ 
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V cos  
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(1.  2) 
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The  requirement  that  the  discrete  potential  energy  function  be  minimized 
with  respect  to  the  undetermined  coefficients  appearing  in  Equation  1.  2 
Results  in  the  generalized  eigenvalue  problem: 


[K  ] | W } + X[C]  {W  } = { o } (1.3) 

Solution  of  Equation  1.  3 yields  an  approximation  to  the  critical  load  and  the 
corresponding  buckled  mode  shapes. 

In  the  computer  program,  at  least  two  solutions  of  Equation  1.  3 are 
performed  per  problem  in  order  to  assess  the  convergence  of  the  numerical 


result.  An  initial  value  of  the  limit  N of  summations  (see  Equation  1.2)  is 
first  selected,  according  to  the  types  of  loading  which  are  to  be  considered. 
Solutions  are  then  obtained  using  N and  (N  + 1)  terms  in  each  series,  and 
the  critical  loads  are  compared.  Further  solutions  are  performed  until  the 
change  in  the  calculated  buckling  loads  is  less  than  0.  250%,  or  until  N 
reaches  a preset  maximum  value. 
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SECTION  II 


PROGRAM  DESCRIPTION 


The  subject  computer  program  (listed  in  Appendices  A-L)  is  coded  in 
CDC  FORTRAN-EXTENDED  language.  The  organization  of  the  program 
and  the  functions  of  individual  subprograms  are  described  in  this  section. 


2.  1 PROGRAM  SEGMENTS 

The  functions  of  each  of  the  individual  sections  of  the  program  are  as 
follows: 


Program  MAIN 

Subroutine  INDATA 
Subroutine  ASSMBL 

Subroutine  BEND1 
Subroutine  CRAMER 
Subroutine  DET 
Subroutine  NROOT 

Subroutine  EIGEN 


initializes  control  parameters,  calls  the 
subprograms,  and  computes  the  execution 
times  of  major  steps  in  the  problem  solution, 
reads  problem  data,  initializes  sizing  para- 
meters and  defines  all  geometric  constants, 
assembles  the  elastic  stiffness  [K]  , and  com- 
putes the  contributions  of  compression  and 
shear  loads  to  the  geometric  stiffness  [G], 
computes  the  geometric  stiffness  terms  due 
to  nonuniform  (edge  bending)  loads, 
solves  a linear  system  of  four  equations  using 
Cramer's  Rule. 

performs  a literal  expansion  to  obtain  the 
determinant  of  a matrix  of  order  four, 
manages  the  solution  of  the  generalized  eigen- 
value problem  [A]  {x}  = \ [B]  {x}  , where  [B] 
is  positive  definite  and  symmetric, 
extracts  the  eigenvalues  and  eigenvectors  of 
a real,  symmetric  matrix. 
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Subroutine  CHECK 


Subroutine  OUTPUT 


Subroutine  ALTER 
Subroutine  TITLE 


locates  a minimum  eigenvalue  and  the 
corresponding  mode  shape,  and  determines 
whether  or  not  the  solution  has  converged, 
prints  the  final  solution  for  critical  loads 
and  mode  shapes.  During  batch  processing, 
intermediate  solutions  and  segment  execution 
times  are  also  printed.  Entry  point  ECHO  is 
called  to  reprint  the  problem  data  for  veri- 
fication purposes. 

accepts  corrections  and  modifications  to 
problem  data  during  interactive  execution, 
prints  a title  heading  at  the  start  of  execution. 


2.  2 PROGRAM  ORGANIZATION 

The  general  structure  and  flow  of  control  within  the  computer  program 
is  shown  in  Figures  4 and  5.  Execution  proceeds  in  exactly  the  same  way  for 
both  batch  and  interactive  processing,  with  the  exceptions  of  a multiple- run 
provision  in  the  batch  mode  and  a data  modification  sequence  for  interactive 
use. 

The  solution  procedure  is  arranged  as  follows.  The  problem  data  is 
read,  and  an  initial  range  on  the  displacement  series  (Equation  1.  2)  is  chosen, 
according  to  the  types  of  loading  which  are  specified  (subroutine  INDATA). 
Subroutines  ASSMBL,  CRAMER,  DET  and  BEND1  are  entered  to  form  the 
elastic  and  geometric  stiffnesses,  and  the  resulting  eigenvalue  problem  is 
solved  in  subroutines  NROOT  and  EIGEN.  The  lowest  eigenvalue  is 
located  »nd  compared  with  previous  solutions  (subroutine  CHECK).  If  the 
required  convergence  test  is  passed,  the  solution  is  printed  (subroutine 
OUTPUT),  and  a new  problem  is  begun.  Otherwise,  the  number  of  terms 
in  the  displacement  approximation  is  increased  by  one,  and  subroutine 
ASSMBL  is  recalled  to  begin  the  next  solution. 
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The  maximum  number  of  terms  which  can  be  considered  in  the 
numerical  analysis  is  limited  by  high-speed  storage  requirements,  since 
no  scratch  files  are  used.  For  interactive  use,  the  program  is  restricted 
to  seven  terms  in  the  displacement  series  (Equation  1.  2),  since  the  resulting 
eigenproblem  is  of  order  hr  = 49.  The  order  of  approximation  can  be 
increased  when  the  program  is  executed  noninteractively,  although  this  is 
usually  not  necessary  for  panels  having  reasonably  small  aspect  ratios. 
Array  dimensioning  for  the  listing  of  the  program  (Appendices  A - L) 
corresponds  to  a maximum  of  N=  7. 


■ 
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SECTION  III 
PROGRAM  USAGE 

The  subject  computer  program  is  operable  from  card  input  (batch 
processing)  or  interactively  at  a teletype  console.  This  section  describes 
the  required  input  data  for  both  modes  of  operation,  and  gives  the  necessary 
information  for  interpretation  of  the  output. 

3.  1 DATA  INPUT  FROM  CARDS 

Punched-card  input  is  arranged  in  free  format,  with  each  item  of  data 
separated  from  the  next  by  a space  or  a comma.  The  data  deck  is  arranged 
as  follows: 

Card  1 : 0 
Card  2 : NRUN 

Data  for  problem  1 

Data  for  problem  2 

• • 

• • 

Data  for  problem  NRUN 

A zero  entered  on  the  first  data  card  indicates  batch  processing.  The 
parameter  NRUN  determines  the  number  of  problems  to  be  considered. 

Data  items  for  a single  problem  are  entered  on  five  cards,  as  indicated 
below: 

Card  (a):  El  E2  NU1  NU2 

Card  (b):  GX  GY 

Card  (c):  T1  T2  TC  A B 

Card  (d):  NX  NXY  NY 

Card  (e):  NXB  NYB 


Some  comments  are  in  order  concerning  the  loading  parameters 
specified  on  cards  (d)  and  (e)  above.  The  values  NX,  NXY,  NY,  NXB,  and 
NYB  are  understood  to  specify  only  the  relative  magnitudes  of  each  of  the 

12 


S 


applied  edgewise  loadings;  the  intensity  of  loading  is  determined  by  the 
eigenvalues  of  the  problem.  For  example,  if 


NX  = 1 . 

NXY  = 2. 

NY  = NXB  = NYB  = 0., 

a critical  load  intensity  is  sought  such  that  N = 2N  . The  values  of 

xy  x 

NXB  and  NYB  refer  to  the  maximum  magnitude  of  an  edgewise  bending  load, 
measured  at  the  corners  of  the  panel  (see  Figure  3).  That  is,  values  of  the 
corresponding  loads  N^»  N ^ are  to  be  interpreted  as  defining  the  edgewise 
resultants 


N = N ( 1 - 2y/b) 
x xb 


N = N . (1  - 2x/a). 

y yt> 


(3.  1) 


The  corresponding  total  bending  moments  are  then  given  by 


M = N . b^/6 
x xb 


M = N L a2/6. 
y vb 


(3.  2) 


The  loading  parameters  specified  on  cards  (d)  and  (e)  may  be  given  either  as 
integer  or  real  values. 


3.  2 DATA  INPUT  FROM  TELETYPE 


The  data  required  for  interactive  processing  is  similar  to  the  card  input 
outlined  above.  However,  the  counter  NRUN  is  not  used,  and  additional 
options  may  be  exercised  for  making  modifications  to  the  problem  data  and 
for  selection  of  output. 

To  initiate  requests  for  input  data,  the  value  "1"  is  entered  (in  free 
format)  as  the  first  item  of  data  following  the  command  to  execute;  for  example. 
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The  program  then  responds  with  the  following  requests  for  input: 


ENTER  MODULI  & POISSONS  RATIOS,  El  E2  NU1  NU2  : 

ENTER  CORE  MODULI,  GX  GY 

ENTER  GEOMETRIC  PARAMETERS,  T1  T2  TC  A B : 

ENTER  LOADING  PROPORTIONS,  NX  NXY  NY  NXB  NYB  : 

After  each  ENTER  request,  the  appropriate  input  data  are  typed  on  the  same 
line,  with  each  item  separated  by  a space  or  comma. 

An  option  is  provided  in  the  computer  program  for  rerunning  a 
previously  defined  problem,  with  one  or  more  input  items  revised.  Following 
the  output  for  the  current  problem,  the  message 

BEGIN  ALTER 
MODIFY  DATA  (Y/N)? 

is  printed.  If  the  next  analysis  to  be  done  differs  from  the  last  in  only  a few 
parameters,  requests  for  revisions  to  the  data  are  initiated  with  the  response 

MODIFY  DATA  (Y/N)?  Y (cr) 

A request  is  made  by  the  program  for  the  first  parameter  to  be  modified,  as 
follows: 

GIVE  * PARAMETER  NAME 

* 

Immediately  following  the  asterisk,  the  user  enters  the  name  of  a program 
variable  whose  value  is  to  be  changed,  and  depresses  the  carriage  return. 


* The  symbol  (^R)  indicates  that  the  carriage  return  key  is  depressed.  All 
commands  and  data  entered  by  the  user  are  underlined. 


Variable  names  recognized  by  the  program  are  given  in  the  list  below. 


El 

E2 

NU1 

NU2 

GX 

GY 

T 1 

T2 

TC 

A 

B 

NX 

NXY 

NY 

NXB 

NYB 

END 

The  program  responds  by  requesting  the  new  value  of  the  data  parameter. 
For  example,  the  process  of  changing  the  modulus  of  the  lower  face  sheet 
of  a sandwich  panel  proceeds  as  follows: 

BEGIN  ALTER 

MODIFY  DATA  (Y/N)?  Y ^R) 

GIVE  ^PARAMETER  NAME 


= 10.  2 E6 


*END 


The  keyword  END  signals  the  end  of  modifications  to  the  data.  Normal 
completion  of  all  data  revisions  is  verified  by  the  message 
ALTER  COMPLETED, 
and  the  solution  to  the  new  problem  is  begun. 

If  no  revisions  are  to  be  made  by  ALTER,  the  appropriate  response  is 

BEGIN  ALTER 
MODIFY  DATA  (Y/N)?  N 

The  next  line  to  be  printed  is  then 

BEGIN  NEW  PROBLEM  (Y/N)? 

The  response  N at  this  point  causes  execution  to  be  terminated. 

A response  of  Y^  initiates  the  input  sequence  for  a new  problem 

beginning  with  the  ENTER  requests  as  previously  described. 
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3.  3 DESCRIPTION  OF  OUTPUT 

Output  from  the  program  for  each  individual  problem  solution  consists 
of: 

i)  input  data  verification 

ii)  results  of  intermediate  calculations 

iii)  problem  solution 

iv)  diagnostic  information. 

The  current  set  of  problem  data  is  reprinted  upon  completion  of  the 
input  procedure  for  each  problem.  During  interactive  use,  the  data  is  also 
printed  after  each  set  of  modifications  in  order  to  provide  a readable 
summary  of  the  current  problem. 

During  batch  processing,  the  results  of  each  solution  iteration  are 
included  in  the  output.  This  information  includes  the  number  of  terms  in 
the  displacement  approximation  (Equation  1.2),  current  values  of  the  critical 
loads,  and  processing  times  for  matrix  assembly  and  eigenvalue  extraction. 
Intermediate  results  are  not  displayed  during  interactive  execution  of  the 
program. 

Final  results  displayed  by  the  program  include  the  critical  load 
magnitudes,  the  eigenvalue  solution,  and  the  mode  shapes  for  the  first 
buckling  mode.  When  executing  the  program  interactively,  the  user  can 
suppress  the  display  of  eigenvalues  and  mode  shapes.  It  should  be  noted 
that,  in  a linear  buckling  analysis,  only  the  relative  magnitudes  of  the  mode 
shape  coefficients  are  meaningful;  therefore,  the  vector  of  coefficients  is 
normalized  to  unit  length  before  printing. 

Diagnostic  messages  output  by  the  computer  program  are  of  two  types. 
In  the  event  that  a solution  using  the  maximum  allowable  number  of  terms  has 
failed  to  converge  within  one  quarter  of  one  percent  (see  Paragraph  1.  2), 
the  message 
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WARNING:  RELATIVE  CONVERGENCE  TOLERANCE 
OF  . 0025  NOT  SATISFIED 
TOLERANCE  OF  ( ) ACHIEVED 

is  printed.  This  situation  will  occur  in  the  analysis  of  panels  having  very 
large  aspect  ratios,  and  subjected  to  edgewise  shear  or  bending  loads.  It 
should  be  remarked  that  the  tolerance  of  . 25%  for  convergence  of  the 
iteration  is  quite  stringent;  in  nearly  all  practical  applications,  convergence 
to  within  one  or  two  percent  can  be  expected  using  a maximum  of  seven 
terms  ( N=  7)  in  Equation  1.  2. 

Messages  may  be  printed  during  execution  of  the  ALTER  sequence, 
to  identify  invalid  parameter  names  (see  Paragraph  3.  2)  entered  by  the 
user.  For  example,  the  following  sequence  occurs  due  to  the  use  of  the 
invalid  name  E3: 

BEGIN  ALTER 

MODIFY  DATA  (Y/N)?  Y @ 

GIVE  * PARAMETER  NAME 
*E3  @ 

INVALID  DATA:  E3 

RETRY  OR  TYPE  END 

* 

Multiple  errors  of  the  above  type  during  a single  pass  through  the  ALTER 
sequence  will  cause  execution  to  be  terminated. 
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SECTION  IV 


1 

EXAMPLE  PROBLEMS 

Typical  analyses  performed  using  the  subject  computer  program  are 
presented  in  this  section.  The  examples  are  intended  to  familiarize  the 
reader  with  the  mechanics  of  input  and  output,  and  to  demonstrate  most 
of  the  important  features  of  the  program.  In  each  of  the  following  examples, 
information  supplied  by  the  user  is  identified  by  underscoring. 

4.1  EXAMPLE  1:  INTERACTIVE  EXECUTION 

This  example  is  an  interactive  session  which  demonstrates  many  of 
the  available  program  features.  Two  separate  problems  are  solved,  each 
involving  various  combinations  of  material  parameters  and  loading  conditions. 

Output  for  the  session  is  shown  in  Figure  6. 

For  the  first  problem,  a 24  in.  by  40  in.  panel  having  0.025  in.  aluminum 
face  sheets  is  analyzed  for  buckling  in  axial  compression.  Following  the 
command  to  execute,  the  problem  data  is  read  in  and  echo-printed  by  the 
program.  A converged  solution,  corresponding  to  three  terms  in  the 
displacement  series,  is  printed  out,  followed  by  the  eigenvalues  and  first 
mode  shape  requested  by  the  user.  It  should  be  noted  that  the  mode  shape 
coefficients  shown  are  not  deflections,  but  modal  participation  factors;  for 
the  present  case,  since  only  W r 0,  the  deflection  of  the  panel  is  of  the  form 

w(x,y)  = A sin  — sin  — u . 

a b 

f 

In  the  ALTER  sequence,  the  face  sheet  thicknesses  are  increased  to  0.  32  in., 
and  the  core  shear  modulus  in  the  direction  of  the  loading  is  made  slightly 
larger.  Once  the  modifications  to  the  problem  are  complete,  the  new  set 

I «’ 

of  data  is  echo-printed,  and  a new  solution  is  performed.  In  this  instance 
the  user  does  not  request  a listing  of  the  eigenvalue  solution  and  buckled 


EH TER  MODULI 


r 


ENTER  MODULI  & POISSONS  RATIOS  * El  EE'  MU l HUE 
c.HT c.F:  LORE  MODULI  - bX  bV 
ENTER  GEOMETRIC  PARAMETER::  - T1  TE  TC  A E 
ENTER:  LOADING  PROPORTIONS-  NX  N' IV  N T N.1E  NVB 


lG.Efcb  l'.i.EEG  .33  . 33 
E 3j  0 0 . £33  'i  Ti . 


i 'JE3  . UC3  #3  E 4 ■ 4 D . 


1 0 ODD 


»( 


FACE  SHEET  PARAMETERS 


FACE  1 : 
F ALE  E • 


E 

. 1 OEGE+OS 
. 1 oe  UE+  ;jy 


COPE  PARAMETERS 


6X 

. esooe+o: 


CORE  : 

PLANFORM  D I MENS  ION S 


H - 
B — 


£4.0 0 0 0 
4 0 . 0 0 0 0 


LOAD  FACTORS 

NX-  1.00  NX 7-  0.00 

MXB-  0.00 


NU 

■ -•  S’  *J  0 *J 

. 33 0 0 0 


, 035 
. «Ju!3 


jfj 

JO 


id  \ 

i330E+05 


N’T—  iJ.DO 
N T B: — 0 . 0 0 


TC 

3 0 0 0 0 


♦♦♦♦♦  SOLUTION  for  M-N-  3 

CRITICAL  LORDS 
NX  - -1168.0004 

NX’T-  0.0000 

N T - 0 . 0 0 0 0 

N X 3 — 0 ■ 0 0 0 0 > v 1 — E ♦’  T ■■■'  B 

N TB  - 0 . 0 0 0 0 * < 1 -£*X  "A ; 


PRINT  EIGENVALUES  AND  MODE  SHAPES  T'T-'N 


*********  EIGENVALUES 

- . 3937EE+04  — .68738E+04  — .64344E+U4  — .33044E+04 
■.44E49E+04  — . 38981 E+ 04  — .33719E+04  — .E6183E+04 

- . 1 1680E+04 

Figure  6.  Computer  Output:  Example  1 
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AT  critical  loads 

•:  W 11'  Id?  . . . 9 W 1 3 9 Wd  1 9 l.'J 2d  ? . . . 9 1'J  33  .* 


1 - 0 0 0 0 • j 0 0 0 0 0 
0 - 0 0 0 0 0 • J 0 0 0 0 
0 . 0 0 • j 0 0 0 0 0 0 ' J 


i.l  . ij  ij  .J  .J  i.l  i.l  ..I  ..I  . j 0 

0 . 0 j 0 ij  o 0 O 0 . j . j 


i.l  . i.1  i.l  ij  i.l  ;_l  .1  i.l  i.l  i.l  i.l 


i.i . i.i  i.i  ,_i  .j  ij  j 0 o .j  .j 


i.l  . '.I  '.I  i.l  ij  'J  i J i.l  'J  ij  i.l 
•j  . j ij  ij  ij  ,j  i j j ij  ij  ij 


BEGIN  ALTER 
NOD  I F V DAT  A *•.  V N 


GIVE  ♦PARAMETER  NAME 
♦ T 1 (CR) 


ALTER  COMPLETED 


FACE  SHEET  PARAMETERS 


7 5 i.i iJ  . CR 


FACE  1 : 
FACE  2: 


.102 0E+ 03 
. 1 <Jd  0E+  iJ  3 


NU  T 

i 3 3 i.l  'J  iJ  ■ ij  -3  d iJ  ij 

1 3-3  'J  'J  iJ  . 03d  0 0 


CORE  PARAMETERS 


CORE  : 


GX 

. c'73  0E+  OS 


. E33  0E+03 


PLRNFORM  HI  MENS  I ONJ. 

A - £4.0000 

E — 4 0 . 0 0 0 0 

LOAD  FACTORS 

NX—  1 .00  NXY-  0.00  NY-  0.00 


MXE-  0.00 


NYE-  0 . O 0 


*****  SOLUTION  FOR  ti-N-  i **** 

Figure  6.  (cont.  ) 
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• * -«wr  r*mTl  t Tj,  t’W  ~ 


*******  CRITICAL  LOAD; 
NX  ■-  -13£1  .3£43 

NXY-  0.  Cl  000 

MY  — 0.0000 

NX 3-  0.0000 

MY 2-  0.0000  * ■ 


1 -£*Y  £ ■' 
i h) 


PRINT  EIGENVALUES  AND  NODE  SHAPES  OiVtOYN 


BEGIN  ALTER 
MODIFY  DATACY-'N.v'N 


BEGIN  NEW  PROBLEM <Y^N)?Y 


ENTER  MODULI  & POISSONS  RATIOS » El  E£  NU1  NU2 
ENTER  CORE  MODULI > GX  GY 
ENTER  GEOMETRIC  PARAMETERS  * T1  T£  TC  A E 
ENTER  LOADING  PROPORTIONS  > NX  NXY  NY  NXB  NYB 


1Q.£3E6  1Q.£3E6  .32  .32 
£ 4 0 00.  £ 0 0 0 0 . 

.063  .063  .375  £0.  £0. 

0 1 0 0 0 ^ 


FACE  SHEET  PARAMETERS 
E 

FACE  l:  . 1 023E+03 

FACE  £ s . 1 0£3E+  03 

CORE  PARAMETERS 

ijX 

CORE  : .£400E+03 

PLANFORM  DIMENSIONS 

A — £0.0000 

B — £0.000  0 

LOAD  FACTORS 

NX-  0.00  NXY-  1.00  NY-  0.00 

NXB-  0.00  NYB-  0.00 

*****  SOLUTION  FOP  M-N-  7 **** 

Figure  6.  (cont.  ) 
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NU 

5£000 
32  0 0 0 


T 

. 063 0 0 
. 063  0 0 


GY 

0 0 0E+ 03 


TC 

.37300 


*******  CRITICAL  LOAD' 


MX  - 
NNY- 
NY  - 
NNB- 
M i 3- 


. 0 0 0 0 

i . 7 ? 1 4 

i . 0 ,j  0 0 

i . 0 0 .j  o ♦ • 

i . o o o o ♦ • 


1 -3  ♦■Y  B 

i : a 


PRINT  EIGENVALUE!:  AND  NODE  SHAPE'S  <Y--N>?Y 


♦♦♦♦♦♦♦♦♦  c.1 

.791 03E+04 
. 1 1431E+03 
. 1 3799E+ 03 
. 1 7443E+  03 


GEN VALUE : ♦♦♦♦♦♦♦♦ 


. 363 1 7E+  04 
. 1 1649E+05 

. 14376E+03 
. 1S433E+03 


1 0403E+0C 
. 1 3 073E+  OC 
1 6 366 E+  oC 
, 3 1 335E+  0C 


. 1 0 3 6 3 E ♦ 0 3 
. 13133E+05 
. 1 663 1 E+  03 
. 33  039E+  03 


S3737E+  03 

.33  i?4r.E+  03 

. 17179E+19 

.S9699E+ 1 9 

34403E+ 1 9 

.43639E+19 

.33136E+19 

.55543E+3 1 

73331E+30 

- . 16733E+30 

— • 1 033 1 E+3  0 

- .37 143E+ 1 9 

33373E+ 1 9 

- . 34393E+19 

- . 14340E+19 

— . 337 46E+ 03 

33737E+ 03 

— . 33  039E+ 03 

-.31335E+05 

- . 1343SE+03 

17443E+03 

— . 1 663 1 E+ 03 

— . 1 636eE+ 03 

— ■ 1 437 6E+ 03 

1 3799E+  03 

— . 1 31S3E+03 

— .13  073E+ 05 

- . 11649E+03 

1 1431E+03 

— . 1 0363E+ 03 

- . 1 04 03E+ 03 

— .36317E+04 

- .791 03E+04 

*********  MODE  shape  ♦♦♦•♦♦♦♦♦♦■ 

AT  CRITICAL  LOADS 

Cl  1 1 ?W13  9 • ■ • j W 1 7 i>W31  ?U33  9 • • . ? W i'' 


87393  0 cl  ? 6 7 

ij 

. 0 0 0 i'i  0 0 0 0 0 0 

.1370574339 

ij  . i.l  iJ  i_l  ij  0 ij  i.l  ij  0 0 

0134394330 

•j 

_ i 1 i i 1 1 i 1 1 1 i 1 i 1 i 1 0 0 

■ ‘J  0581431 5 3 

i.l  . i.l  1 1 i I 1 1 1 1 i 1 i 1 t t i 1 PI 

4334031377 

a ij  ij  i 1 i 1 1 1 i 1 i 1 ij  i.l  0 

- . 0309763331 

0 . 0 o 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 

0 066973347 

:"i 

. 0 0 0 0 0 0 0 i'i  0 0 

. 1 3333  0 0 1 99 

i.l  . ij  iJ  i.l  iJ  iJ  i 1 iJ  i.l  i 1 iJ 

1413367397 

:“i 

. 0 0 0 ij  0 ij  i.l  i.l  ij  0 

-.01 0 1 696369 

0 . 0 0 0 0 0 0 0 0 0 0 

0 J 1 07234  0 5 

0 

■ i J 0 i 1 0 i.l  i_l  • 1 1 1 i.l  i.l 

- . 0355943697 

i.l  . 0 i.l  iJ  i.l  • 1 i 1 i 1 i_l  1 1 i 1 

0373419133 

ij 

_ i 1 fl  i 1 i 1 i 1 i 1 i 1 i 1 0 • J 

. 0069832008 

ij  . i.l  i.l  i.l  ij  ij  i.l  ij  0 ij  0 

0 1 1 433697  3 

•J 

■ 'J  iJ  iJ  'J  ij  'J  0 0 0 U 

— - i.l  1 i.l  d 6 46  i.l  8 ij 

i J , iJ  iJ  0 i 1 i 1 0 0 i J 1 1 i 1 

0 U 66047 0 y 3 

:“i 

. 1 1 1 1 i 1 1 1 1 1 1 1 0 ; M 0 

— . 0 Onfc'CC'  !o4  1 

ij  . i.l  i 1 ij  i 1 i 1 i 1 i 1 i 1 i 1 i 1 

0 o52  ij  7 9 iil  08 

•J 

. 0 0 i'i  i'i  0 0 i'i  i'i  0 0 

• i.l  0 6 J J 3 il!  6 1 

0 . 0 o O i l i 1 0 i l ; l i lil 

0 0 8197  o 3 0 

ij 

. 1 1, Iil. 1.  i.l.  I(|.  Iil 

. 0 03363  i ■'  33  1 

0 . 0 0 • 1 0 1 1 1 1 ; i; i , i , I 

O 0 1 63 28 048 

ij 

. 0 0 i’i  0 i'i  o . 1 i*l  0 0 

— . ij  Oc’8  75  y 1 45 

•J  . i.l  i.l  0 1 1 i 1 i 1 i 1 i 1 0 i 1 

0 * J iil  4 c!  9 4 065 

BEGIN  ALTER  . 

MODIFY'  DflTAvY  'N>  7 Y_  ^ 

GIVE  ♦PARAMETER  NAME 
♦NXY 


NXY  - 0. 

♦ NY  6r) 


Figure  6>  (cont.  ) 
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r 


NY 

♦NXB 
NX  3 - 1 

♦END 


filter  completed 


FACE  SHEET  PARAMETERS 


FACE  1 : 
FACE  2: 


. 1023E+03 
■ 1 US3E+ US 


CORE  PARAMETERS 


GX 

2400E+05 


NU 

i S c!  0 U U 
, 32  0 0 ft 


, U 3 U J 
, 063 0 0 


GY 

O00E+03 


CORE  : 

PLANFORM  DIMENSIONS 

A — 2 0 • 0 0 0 0 

E - 2 U - 0 U 0 0 

LOAD  FACTORS. 

NX—  0 • 0 0 NXY—  0.00  NY—  3.00 

N X B — 1 . 0 0 N rE-  0 • 0 0 

♦♦♦♦♦  SOLUTION  FOR  M-N-  6 ♦♦♦♦ 

♦♦♦♦♦♦♦  CRITICAL  LOADS  ♦♦♦♦♦♦♦ 

NX  - 0 . 0 0 0 0 

NXY-  0.0000 

NY  ^ -3149.7133 

NXB-  -1716.3713  ♦ 1-2*Y-'B> 

NYB-  0.0000  ♦ *■.  1 — S ♦ X y A 1 


TC 

1 3 7 3 u U 


PRINT  EIGENVALUES  AND  MODE  SHAPES  <Y-N>?N 


Figure  6.  (cont.  ) 
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BEGIN  HLTER 
MODIFY  DhTFi.:  Y..-N  :.  tm_ 


BEGIN  NEW  PROBLEM 


STOP 

11.928  CP  SECONDS  EXECUTION  TIME 


Figure  6.  (concluded) 


5 


■ 


The  second  problem  in  the  session  consists  of  a square  panel  with 

identical  face  sheets  and  orthotropic  core  (G^/G^  = 1.  20),  analyzed  for 

buckling  under  two  different  loading  conditions.  First,  a shear  loading  is 

considered,  for  which  a seven-term  solution  is  required  to  obtain  a converged 

result  (it  can  be  seen  that  the  shear-buckling  mode  shape  of  the  panel  is 

quite  complicated;  thus,  more  terms  are  needed  in  the  displacement 

approximation).  The  ALTER  sequence  is  next  used  to  modify  the  data  for 

the  second  loading  condition.  A critical  load  magnitude  is  then  found  such 

that  N = 3N  , , where  N . defines  an  edgewise  bending  load  as  defined  in 
y xb  xb 

Equations  3.  1 and  3.  2. 


The  session  is  terminated  when  the  user  bypasses  the  ALTER  sequence, 
and  does  not  elect  the  option  of  beginning  a new  problem.  Execution  time 
for  the  session  is  displayed  following  the  normal  exit  from  the  program.  It 
should  be  noted  that  the  solution  for  a compressive  load  only  is  started  with 
two  terms  in  the  displacement  series,  and  for  bending  and  shear  loads  with 
five  terms  in  the  series.  Thus,  the  execution  time  listed  for  this  session 
(11.928  sec.  ) is  the  time  required  for  nine  separate  eigenvalue  solutions, 
even  though  only  four  separate  problems  have  been  solved. 


4.  2 EXAMPLE  2:  BATCH  MODE  EXECUTION 

Two  problems  are  solved  in  this  example  to  illustrate  the  use  of  the 
program  in  the  batch  execution  mode,  with  data  input  from  cards.  Output 
for  the  problems  is  reproduced  in  Figure  7.  Input  data  required  for 
running  both  problems  in  a single  submission  is  as  follows: 

Card  1:  0 


Card  2: 


2 


(Flag  for  batch  mode) 
(Number  of  problems) 


c.occe  •<i-2»x/a> 


*********  iI9iNVAL'J£S  ******** 

- . 2ST?3r+  - ; ir« gr+  rr  . ra?c sz* <m  -.ms^au 

-.G&GGlc+CG  •*.45132£>'.'4  -.29892£+C4  1 92 9 1E+ 04 


inrrr  shsps"  »*»»»*»»» 
AT  CRITICAL  LOAOS 


( Mil, W12, .. . ,W13 , W21.W22  . .. . , «33> 


1.00030; j 3 0 c 

3. 33-333333  33 

o.oooooooooc 


o.ooocccocoo 

C • 3 33T333C  j 


0. COCQJCJGOO 
3 . 03  u u 3 u 0 3 33 


o.  a 1.000C  oooo 

3.C3  J3C33033 


Figure  7.  Computer  Output:  Example  2 


mm ipm  y ,'%n 


Mnl 


FACE  SHSHT  PARAMETERS 


PACE 

li  .10001+08 

tra 

. 3 0 C 0 C 

T" 

• 05  00  0 

FACE 

21  • 1 0001  + C 8 

. 3 300  3 

■ C5  000 

CORE 

“ARAMcTERS 

CORE 

‘ ' “5X 

I .25031  + C 4 

~5  y 

• 25  0 3E  *04 

re 

. 25000 

RLANFCR*  DIMENSIONS 

As  1 5.  0 0 : 1 

L 040  FACTORS 


NX  = 

.53  NX  Y = 

3.30  NY=  ,53 

NXB=" 

T.Tr- 

ftyoI'I,  c ; 

TERMS 

NX 

NX9  NX  Y 

NY 

NY9 

ASSEM(SEC) 

EIGEN  ( SECt 

5 

6 

-143.558 

-143.769 

-287.115  -861.348 

-281. 538  -844.614 

-143. 558 
-140.769 

-574.232 
-563#  076 

.355 

.095 

1.530 

4.520 

7 

-139.53  3 

-2797275“  -5  37.  327 

* 13-* 6 ;5 

-P3S, 55Z 

.15  9 

:i,4Co 

»*  EARNING  ! 

of  . 30  25 
»•  TOLERAMC 

R1LAT  I‘/2"  CON  9130i  NCI 
NOT  SATIScIiO 
1 OF”  .0081  ACHIEVED 

TOLER ANC 

r ■ 

+wwww  5Q LTJTTON  FOR  M=N=  7 **** 

******  *~CRTTICfli:  TiCTfifTS-  ******* 

NX  = -139.6379 

NXY  = 

NY  = 

-837. 8274 

-139.6379 

FTTHi 
NY  9 = 

-2  79.2  75  3 '"^Tl 
-558. 5516  * ( 1 

- Z^Y/TJ 
-2*X/ A) 

IGENVAlUtS  ******** 

• 39b 11 w*  33 
•58 136E+G  3 

• 4wZ3bc  *U3 

• 587 15E*C  3 

• 4 6 j 

.71136EK  3 

• 47Q3  3i * 0 3 
. 72757E+03 

7 79  390 1 + 0 3 

. 13640  E+  04 

. <3296‘+c  *0  3 
. 1395 1£*C  4 

•923b3E+C3 
• 1 7 94  9t  *0  4 

. 1UU  331  + 07 
• 1818  9E+  04 

•42521E+34  . 454351+04“ “ rii‘677r+C5  -.9725  81+04 


-.425981+04  -.344011+04  -.263951+94  -.245561+34 
-. 2 3 551 2+  34  -.21  I55T+OV+X1 3I54E‘+Ol.  -.16J941+04 
-.141631+0  4 -•  114  26»+  34  -.  10  0736  + , 14  -,92326c+03 

-.851771+03  -.  681781  + C~3~  - . 565061  + 4 3 -.622591+03 

-.633331+  03  —.  528351  + 03  -.524241+C  3 -.  463631  + 03 

+9 4 58+81+ 03~~=Y  4020  32*23 ~ 398921  + CT  - . 3772 41+27 

-.37313E+03  -.311721+03  -.308951+03  -.283931+03 
+.279281+03 


' + + wool  SHAP1"**» 
AT  CRITICAL  LOAOS 


(W11.M12.... 

.H17.W21.H22 

W77) 

-.3122565  41J 

• **83833221 

. 0234377730 

3504191233 

.1389685001 

. 0233234855 

-. 024C 176481 

.0392184216 

.1203941242 

-. 4678519343 

. 2304119503 

. U34325660 

— . 15  31831671 

• 0 -*6 - 8o955  a 

• j 5221  74 u 11 

-.1583353  382 

.10  178  33486 

. 1633781658 

-.  231260500  4 

. 0 789717  120 

.0197829947 

-.  Ol,  i216552o 

. UU33624638 

. 04^1633053 

-.1  J04224132 

. 0 37210520  7 

. 0451  371127 

-.036+715153 

•0023476607 

+ .807X078828 — 

-• CIO  9 879394 

• c zt'i'i  4 J 3 it 

.00 12483311 

-.  0 1974,3:318 

. 0048219853 

.3004122317 

-.0003442618 

. 0C874i*77u 5 

-.0117641349 

-.0549153268 

.0117571273 

-. CO  39650  4 7 3 

. 0006  310  114 

-. 0121763316 

-.0020389  207 
-«0i,JS36O2C6 

. OC  44697813“  ‘ 

. 0025877969 

“ -.0035189119 

Figure  7.  (concluded) 
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Card 

3: 

10.  2E6  10.  2E6 

.30  .30 

Card 

4: 

23500.  23500. 

Card 

5: 

.023  .023  . 375 

23.5  23.5 

(Data 

Card 

6: 

1 0 0 

Card 

7: 

0 0 

Card 

8: 

10.  0E6  10.  0E6 

.30  .30 

Card 

9: 

2500.  2500. 

Card 

10: 

.050  .050  .250 

15.  30. 

(Data 

Card 

11: 

.50  3.0  .50 

Card 

12: 

1.0  2.0 

Card 

13: 

7/8/9  END  OF  RECORD. 

The  first  case  considered  is  that  of  a square  panel  with  an  isotropic 
core,  analyzed  for  buckling  under  a single  compression  load.  It  can  be  seen 
from  Figure  7 that  the  output  is  identical  to  that  obtained  in  an  interactive 
session,  with  additional  information  printed  concerning  intermediate 
calculations  and  execution  times.  For  the  present  problem,  the  solution  is 
begun  with  two  terms  in  the  displacement  series,  and  converges  immediately 
due  to  the  simplicity  of  the  buckling  mode. 


The  second  problem  considers  a panel  with  a 2:1  aspect  ratio,  subjected 
to  a complicated  system  of  loading  which  includes  compression,  bending 
and  shear  forces.  All  input  data  is  echo-printed,  followed  by  the  results  of 
intermediate  calculations.  Due  to  the  complexity  of  the  loading,  the  solution 
is  begun  using  five  terms  in  the  displacement  series,  and  proceeds  through 
the  maximum  number  of  terms,  (seven).  Since  the  relative  change  in  the 
last  two  solutions  is  greater  than  the  preset  tolerance  of  0.  25%,  a diagnostic 
is  printed,  followed  by  the  current  solution.  In  this  case,  the  diagnostic 
indicates  that,  although  the  solution  has  not  fully  converged,  the  relative 
change  in  the  result  is  less  than  one  percent;  hence,  it  appears  that  the 
buckling  solution  is  still  quite  accurate. 
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4.3  EXAMPLE  3:  A DESIGN  PROBLEM 


This  example  illustrates  the  use  of  the  program  for  a simple  design 
application.  Consider  the  design  of  a simply- supported  sandwich  panel  which, 
due  to  design  constraints  and  material  availability,  has  the  following 


properties: 

Face  sheet  material 

Core  shear  moduli 

Core  thickness 

Planform  dimensions 
Maximum  thickness 

Minimum  thickness  of  faces 


- Aluminum  5052-H32 

E = 10.  2 x 10&  psi,  v = 0.  33 

- G = G = 21300.  psi 

x y 

- t = 0. 5 in. 

c 

- a = 30.  in.  b = 24.  in. 

- t,  + t_  + t s 0. 65  in. 

12c 

' ti  = tz  ^ 016  in- 


The  face  sheet  thickness,  t,  is  to  be  selected  from  standard  gage  sizes 
(.016,  .020,  .025,  .032,  .040,  .050,  . 063)  so  as  to  satisfy  stability 
constraints  for  the  following  loading  conditions: 

(i)  N = N = -750.  lb/in 

x y 

(ii)  N = 600.  lb.  /in.  ; N , - 300.  lb/in. 

xy  xb 

Output  for  the  present  problem  is  shown  in  Figure  8.  A solution  is 
performed  first  for  the  compressive  loading  case,  using  the  minimum  gage 
size  of  .016  for  each  face.  Load  condition  (ii)  is  next  checked,  to  determine 
which  loading  case  is  more  critical  in  determining  face  sheet  thicknesses. 
Since  the  panel  is  unstable  only  for  load  case  (i),  the  facing  thickness  is 
increased  to  . 020  in.,  and  both  load  conditions  are  reanalyzed.  The  critical 
loads  are  then  found  to  be 


(i)  N = N = -809.  lb.  /in. 
x y 


(ii)  N = 3249.  lb. /in.  ; N ,=  1624  lb.  /in. 
xy  xb 


and  the  proper  face  sheet  thickness  is  therefore  t = .020  in. 
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ENTER  MODULI  & POISSONS  RATIOS?  El  E2  NU1  NU2 


1 0 . LE6  10.2E6 


ENTER  CORE  MODULI?  GX  GV 

ENTER  GEOMETRIC  PARAMETERS  ? T1  T£  TC  A B 

ENTER  LOADING  PROPORTIONS?  NX  NXV  NY  NXB 


FACE  SHEET  PARAMETERS 


E 

FACE  1 • . 1 02OE+03 

FACE  2 : . 1020E+03 

CORE  PARAMETERS 


NU 

. 33 0 0 0 
. 33 0 0 0 


GX  GV 

CORE  : .2130E+05  .2130E+05 

PLANFORM  DIMENSIONS 

H — 3 0 . 0 0 0 0 

B — 314 . 0 0 0 0 

LOAD  FACTORS 

NX-  1 .00  NXV—  0.00  NY—  1 .00 
NXB—  0.00  H'i'B-  0.00 


SOLUTION  FOR  M-N-  3 ++++ 

♦♦♦♦♦♦♦  CRITICAL  LOADS  ♦♦♦♦♦♦♦ 

NX  — —646.3333 

N.i  r — 0 ■ 000  0 

NV  - -646.3333 

NXB—  o . o 0 o 0 * v 1 — 2 B !' 

N't'B-  0 . 0 0 0 0 ♦ •:  1 -2*X-A  > 


PRINT  EIGENVALUES  AND  MDDE  SHAPE! 


BEGIN  ALTER 

MODIFY  DATA' V N>?V_  ^ 

Figure  8.  Computer  Output: 


■ 3 1 3 0 0 . 3 1 0 0 . 


• . 0 1 6 .016  .3 
NYB  : 1_  0_  1 0 Q (SR) 


T 

ij  1 'r>  ij  ij 
ij  1 A ij  ij 


TC 

.50000 


(V-  N)?N_  @ 


Example  3 


GIVE  ♦PARAMETER  NAME 
♦MX  PE) 


NXY  - 1 


♦MX  3 


NXB  - .r 


♦END 


ALTER  COMPLETED 


FACE  SHEET  PARAMETERS 


FACE  1 : 
F ACE  d • 


. 1 iJE  0E+  Liy 


CORE  PARAMETERS 


T 

• ij  1 G>  ij  ij 


1 '.I  C ijE  ♦ ijy  • 33  »j  U J . i.l  1 r.  i_l  ij 


CORE  : 


GX 

.cl  3 'JE+  '.'3 


. d 1 3 uE+ >j3 


TC 

3 0 0 0 0 


PLANFORM  D I MENS I ONS 

A - 3 0 . 0 ij  0 £i 

3 — c4  . ij  0 0 ij 

LOAD  FACTORS 

NX-  0 • U i j NXY-  1 .00  M l'—  ij.ijij 


NX 3-  .30 


NY  3—  iJ.iJiJ 


*****  SOLUTION  FOP  M-N-  6 ♦ ♦♦♦ 


Figure  8.  (cont.  ) 


♦♦♦♦♦♦♦  CRITICAL  LOAD  I-  ******* 


NX  - 

0 . 0 0 0 0 

NXY- 

— b!6Gb  . 477 0 

MY  - 

l_l  a i_l  i_l  iJ  tl 

NXB- 

— 1 1 j “I ! -1  m q *!  3 5 

♦ <.  1 -c>Y 

•••■'  B 

NYB- 

0 . • j 0 0 0 

♦ 1—c+X 

.-■R 

P PINT  b I b c. M V A L J E . A H D M DDE  _HhPE S ' r N > 7 N 


BEGIN  ALTER 

MODIFY  DATA<Y/N>?Y  g 

GIVE  ♦PARAMETER  NAME 

♦11  ( 

Tl 

— . Ob!  0 (C 1^) 

♦ T£  { 

g 

T£ 

= .020  g 

♦END  ( 

g 

ALTER  COMPLETED 


**•> 


FACE  -SHEET  PARAMETERS 

E NU  T 

F ACE  1 i . 1 Ob!  OE  + US  . 1 i 0 iJ  0 .Ob!  0 U 0 

FACE  £ « .1  Ob!  UE+  08  . 33  0 0 0 « ObOOO 

CORE  PARAMETERS 

GX  GV  TC 

CORE  i . 2130E+O3  .213OE+05  .30000 

PLAMFORM  D I MENS  I ONS 

A — 3 0 . 0 0 0 0 

B — £4.0 0 0 0 


Figure  8.  (cont.  ) 

X 


M 


FACE  IHEET  PARAMETERS 


FACE  1 : 
FACE  £: 


. 1 Oc  OE  +•  03 
. 1 OE  0E+  03 


MU  T 

■ J -•  0 ij  ij  • i.l  zl  0 ■ J ij 

■ -•  » J • j . i.i ij  i.l  i.l 


CORE  PARAMETERS 


CORE  : 


bX 

.£13  OE+  OS 


GV 

: 1 3 0E+ 05 


PLRNFQRM  DIMENSIONS 

H — ! Zf  ij  , ij  U »J  ij 

? — c!4  a i.l  i.l  »J  ij 

LORD  FRCT0R3 

NX—  1 . 00  NX' f'—  0 . 0 0 NY—  1 .00 


NXB—  0.00 


NYB—  0.00 


▼ SOLUTION  FOR  M-N-  3 ■**++ 

CRITICRL  LORDS  ♦♦♦♦♦♦♦ 
NX  - -309.1307 

HXY  — 0 « 00  oO 

NY  - -309.1307 

NXB—  0 • 0 0 0 0 ♦ 1 — £>’i' B .> 

NVB—  0.0000  ■*  1 — £ ♦ X R 


PRINT  EIGENVALUES  RND  NODE  SHAPES  CY-TO  ?N 


BEGIN  ALTER 
M 0 D I F ' i'  DAT R ( ' i •••' N > T" N 


BEGIN  NEW  PROBLEM < Y.- TO ‘N 


STOP 

1 £ . 0 1 1 CP  SECONDS  EXECUTION  TIME 


Figure  8.  (concluded) 
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SECTION  V 

SUMMARY  AND  RECOMMENDATIONS 


A computer  program  for  the  analysis  of  stability  of  flat,  rectangular 
sandwich  panels  has  been  presented.  The  program  is  capable  of  considering 
very  general  systems  of  applied  loading,  and  takes  into  account  the  effects 
of  orthotropic  core  materials  and  unbalanced  face  sheets. 

The  subject  computer  program  is  capable  of  obtaining  solutions  with 
a relatively  small  amount  of  computation,  and  requires  very  little  computer 
storage  (the  compiled  program  occupies  11460  decimal  words  of  memory). 

As  such,  it  is  ideally  suited  for  use  in  optimization  or  automated  design 
systems. 

A number  of  extensions  of  the  present  analysis  are  possible,  and 
should  be  investigated.  Within  the  confines  of  a linear  analysis,  generalizations 
of  the  present  code  to  include  layered  composite  faces,  elastic  edge  rein- 
forcements and  more  general  panel  shapes  would  represent  a valuable 
capability.  The  existing  code  can  be  extended  to  consider  the  case  of 
multicore  constructions  and  panels  which  undergo  normal  deformations, 
although  both  of  these  modifications  require  the  solution  of  larger  numbers 
of  degrees  of  freedom.  The  analysis  of  other  than  simply  supported  panels 
can  be  achieved  by  altering  the  assumed -mode  approximation,  although 
computational  effort  will  be  increased  for  other  cases.  * 

In  conclusion,  a useful  and  reliable  tool  for  the  analysis  of  sandwich 
structural  components  has  been  described  herein.  Furthermore,  the  subject 
computer  code  provides  a suitable  starting  place  for  the  consideration  of  many 
more  general  types  of  sandwich  construction. 


* 
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SUBROUTINE  -I'ns’i 

^gAL  LX,LX3,LY,LYS,LXY,NX,NXB,NY,NY3 

cerntON  /-Tft  Trrr  r?r,  rt^-rr xmrrZTrx  tntrtt  m-r^r  xtzt  I wrzn 

com^ion/:qr3/gx,gy  , tc,a  , b 

COMMON /O NT L/TNTEp ,PR£V ,VAL  . rSTART, K,NAX,  NRUN,  ICONV  , TIHE1  , Tiweg 

C0MM0N/LCAD/LX,LX6.LY,LYa,LXY,NX,NX3,NY,NYg 

ttWTCA-t-TinS^ 

IF < I NT  ER  > GG  TO  2 

TF-rrsrATT  rtrzrrt 

1 P.5A0  (5  ,*  ) NRUN 

— ? ItWVrfl 

MAX=7 

Pft£\F=0 ~i 

LIMIT=13 

TF  tMRtmrG  T , tT  HIT  rS  TOP 

IF ( I NT  S 3 ) WRIT  I (6  , 100) 

RgAf>Y5  , -)  £ (1>  , g(  2>  ,XNU<1 ) ,XNtlfgi 

IF  { INT  FR  > WRIT  E (6 ,125 ) 

M40  XrG  Y- 

IF (INTER)WRIT  £ (6,150) 

RgAfH5  ,'*'T~T  

IF (INTER)WRITE (6,175) 

RE  AO  (5  , ‘ HfX-ft-X-Y , N Y',  V Xe^ltYfl 

ENTRY  RESET 

-t-X  B=-2*^  X 3 

LY  0=2*  NY  3 

t*-=NX— N-X^ 

LY  = NY -NY  3 



IF (LXY.Nc. 0)K=? 

I fHtXfl  W=»-9-.E>RvL-YB  <=5 

THICK=T(2) +TC  + T < 1 ) 

RIGN-(  1 )■=— -1-, 

SIGN (2 ) = 1 . 

— « e-3-  -f^t-r  2 

XLMOA ( H =1,-XNU( I » *XNU  (I) 

3 RRT  A ( I >=-T  fl-)— THICK /2.- 

TSTART=1 3 

RETURN 

130  PO-HAT  (/2X  ,A3HFNTcR  MODULI  POISSONS  RATIOS,  El  E2  NU1  NU2  t ) 

— 125 FCRM-A^T-t— EX-^-YgN-T^gR— eeRE'  -HOOWrl-, 6-X GY — * 

1 1 3 X , 2H  ; ) 

15  0—  FD°M  A'Y  f--2X-,4«3-HE,'H5R—GEOMETRIG— PAR  A METERS, — T-l— T-2— T€— A— 8 , 

1 6X,2HS  ) 

1T5-  eOPM4-T  f-  >X  , A3H-ENTER-tOACMNG  P^ROPORT-IONS  , NX  NXY--NY-NX0-  NYR — 1-1- 
END 


APPENDIX  C 
SUBROUTINE  ASSMBL 


*: 


s 

I 

* 


\ 


1 
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SUGROOTNE  ASSMGL 

Pr.A L L X , L X 3 , L Y , L Y B , l X Y ,NX,NXB*NY»NY3 

— otr*n  C'r/M-a  t l at  1 2 r , t 1 2 r v an  q-  (-2rr  al  rms- 1 ?>-,  be  t a ( z a , s i g n 1 2 r 

COMMON  /:  OP  E/G  X ,G  Y , TC  ,fl  , B 

^a-^Ton-/^yTirArnj*2^vPT2^ratr7~rrT-ap^pc,^aA , t PtrN-rt^onv , t i we  i ,-t  ingg- 

COMMON/SCLN/A A (24C1) ,B8 (2401) ,EIG<49> , TIGV (2401)  * NOTE 
C 0 CN/t.  C A e/t  A, L X 9 , L f-f  t A3"»  LAY-,  NX-,  NXB , W Y , N Y 9 - - 
•DIMENSION  AC  ( 4 »-»  ) » BC  ( 4 ) , XC ( 4 ) 

L OG  IOAtr  I NTS  9" 

DATA  =1/3.1415325535393/ 

~X2^K*X 

<4=K2*K2 

— 1 DO'  1-1=1 , <9 

A A (I)  = 0. 


DO  110  N=1,K 

">0—  1'tO'  -t-f,* 

JJ=K*(M-l) +N 

JJ-tA2*“f  i J-D-t- J-J 

GO  3 1=1,4 

GO  2 J»±rte 

2 AC  ( I , J ) = 0 . 

-3 

0 M = M*P 1/ A 

<>=«•»=  T/B — 

GO  50  1=1,2 

01=5  M I*-? ETA  M 1 AXLMOAfT) 

C2=E (I)*T(I)*DETA(I)*SIGN(I)/XLMOA(I) 

0-T=E-tI  r* T TlT-AXtMtGft-TM 

C4  = C3*TC I) *SE  T A(  I ) 

05=0  3*-TH>-«-ST9MI> 

AC (1 , 1 ) = AC  ( 1 » 1)+ ,5*GX*TC+C1*  CCM*CM  + ,5* (l.-XNU(I ) ) *CN*CN> 

A C ( 1 >2- ) = AC  M,  2 H-.  5*C1*  ft.-t-XNUM T1*GM*CN 

AC(1.3)=AC(1,3)+C2*(CM*CM+.5*(1.-XNU(I))*CN*CN) 

AC-tl-r^  Y^A  C- ( Tr— Y-^irT-’-e  2*€tt*€ TTr^ Yt+thTT T ) 

AC (2,2 ) = AC (2,2)*.  5 *G  Y* TC +C1 M CN* CN  + . 5 * (l.-XNU(I))*CM*CM> 

AC  f 2 , ^ >=AC  (2,41+0  ?*f  GH*CH+.  9+M-.  -XNUH  H *Cu*CMD 

AC (3,3)  = AC  <3,  3>+C3*(CM*CM+. 5* (l.-XNU(I) ) *CN*CN) 

A C < 3 , * 1 = AC-  T3^  >►-.  5 *C  3 ^ U . N O < T >-  > *0  M*0  H 

AC  (4,4 ) = AC (4, h ) + C3* (CN*CN+. 5* (l.-XNU( I )) *CM*CM) 



DC (2)=9D 12) +. 5*0“* CM* ( CN *CN+CM*CM ) +GY*TC*CN/2. 

90  (3  V = 90  < 3 >■+.  = *6  = *CM*  (CM*BH-+eH*C!4) 

5 0 BC  (4)  = B0  (4)  + . 3*C5*CN* (CN»CN+CM*CM) 

GO  95~  1 = 1,4 — 

55  BC ( I ) = - 3C ( I ) 

AOt  2rt-)==+e-tl-s-2J) 

AC (2 , 3 )=AC (1. 4) 

A C C3^1  >-*C  (1,-3  ) 

AC(3,2)  = AC (2,  3) 

AC(4,1  >=AC  (1,  •>> 

AC  (4.2)  = AC (2,  4) 

ACT4T3  )=ae(-3-,-4T 

305=3. 

GO  ?D  1=1,2 

C3  = E ( I ) * T ( D/XLYQA  (I ) 

70  3CC=BC5+(C3*''  ( I ) * T (I)/3.  ) * ( CM*CM+CN*CN  ) *-  ( C“»CM+CN*CN  ) 
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BCr=0C5+CM*CM*GX  *TC+CN*CN*GY*TC 
GALL  CcAMER(AC,BC,XC>  '' 

A A (J  J)  =0C5 
TO  133  1=1,4 

130  AA(JJ>=AA(JJ)-PC(I)*XC<I> 

110  CONTINU* 

00  15  0 *1=1,  K 
00  153  M=1,K 
00  150  I R=  1 , K 

00  150  IS=1,< 

CM=H*PI/A 

C N=N* P 1/ 3 
IRf1W=(M-l)  *K+N 

1 COL  = ( IP- 1 ) *<  + IS 
JJ=(ICCL-1 »*<2+IR0W 

IF  (I  ROW.  EQ.ICODGC  TO  125 

ICHK=(«+IRI/2 

XCHK  = . 5*  F L CAT  (M+IF ) 

IF (AOS  (FL0 AT (ICHK ) -XCHK ) ,LT, • 100)GO  TO  150 
TCHK= < N+IS ) /2 
XCHK=.  5*FLCAT  CN+IS) 

IF  (AOS  (rL0 AT  ( ICHO-XCHK)  ,LT . , 100  ) GO  TO  150 

OB1JJ) =-32.*LXY*M*N*IR*IS/ (A*B* < IR*IR-M*H ) * (N*N- IS*IS ) > 

GO  TO  150 

125  OB<  JJ)  =-CH*Cwl*LX-CN*CN*LY 

153  CONTINUE 
RETURN 
ENP 


APPENDIX  D 


SUBROUTINE  CRAMER 


I 


?» 


L. 
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SUBROUTINE  CRATER  < A , O , X ) 
DIMENSION  M'4,4|,B(u),X(4),T(4) 
CALL  OFT  ( A ,OENOM  ) 

00  BOO  J=1 ,4 
on  450  1=1,4 
T (T)  =fl (I , J> 

450  4<I,J)=0(I> 

CALL  OFT(A,V(J)) 

X ( J) =X ( J) /OENOM 
00  475  1=1,4 
475  A(I,J)=MI) 

500  CONTINUE 
R E T U R N 
FNP 
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APPENDIX  E 


SUBROUTINE  0;.T(A,P) 
niM"N5liN  a<h,4) 

ni=(A(3,3)+A<4,4)-A(4,3)+A<3»4))+<A<l,l>+A<2,2)-A<l,2>+A(2tl>) 
02=(A(3,2)*A(^,4)-A(4,2)*A(3,4)>4(A(1,3)*A(2,1)-A(1,1)*A(2,3)) 
0 3=(A(3,U*A(+,4l-A(4,H*A(3,4n*(A(l,2)*A(2,3)-A(l,3)*A(2,2n 
04=<A(3,2)+AU,3)-A(4,2)+A(3,3nMA(l,i>*A(2,4)-A<l,4>*A(2,i)> 
nc=(A(3,l)*A(^,3)-A(4,l)»A(3,3))*(A(l,4)*A(2,2)-A(l,2>*A(2,4)) 
06=<A(3,l>*A<4,2»-A<3,2)+A(4,l>)+<A(lt3>*A(2,4)-A(l,4>*A<2,3)) 
0=01+02+03*04 + OS +06 
Tu=>N 

END 


1 

* 

it 


ZA 


SU330'JT:nE  N3OOTCM,A  ,3, XL, XI 


uctaa  ii  ■ ■■  leanc  itjiiii 


pu3Dosr — 

COMPUTE  rI GENV ALUF S AND  EIGENVECTORS  OF  A REAL  NONSYMMETPIC 

Ml  TR  It-  OF  THE"  FORM  S"- 1 NVfRSC"- TIKES-  tc. T HI  S"  SU8EOUT  INC  "IS~~ 

NORMALLY  CALLEO  BY  SUBROUTINE  CA NOR  IN  PERFORMING  A 


- -trsA-o  

CALL  NROOT  (M,A,B,XL,X) 


DESCRIPTION  OF  PARAMETERS 


A - INPUT  MATRIX  (M  X M>. 

C INPUT  M-frTRI  X — fM— X~H>t 

XL  - OUTPUT  VECTOR  OF  LENGTH  M CONTAINING  EIGENVALUES  OF 

NVE  PSE-  TIMES  "A-i — 

X - OUTPUT  MATRIX  (M  X M)  CONTAINING  EIGENVECTORS  COLUMN- 


IP  A DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  OESIRED,  THE 
C I N-COtUMM-i-  GHOULO  ae-PEMO-VeO-FRQ+F-T-HE-  OOUBLE-  PRECISION  — 
STATEMENT  WHICH  FOLLOWS. 


DOUBLE  3 RE  CIS  I ON  A , B , X L , X , SUM V 


THE  S must  ALSO  BE  REMOVED  FROM  OOUBLE  PRECISION  STATEMENTS 
APPEARING  IN  OTHEO-ROUTINES  US  CO  “IN  CONUUNC T I ON-W1TH-THI S 
30U  T IN  E . 


THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 

-0 ON T A I N- DO NBtr£-PR =e IS-I ON-POPTR* N-PUNC  T-fONS^ S ORE — IN— ST-ATEME NT  S 

110  AND  17"  MUST  BE  CHANGED  TO  OSQRT.  ABS  IN  STATEMENT  110 
MUST  BE- CHANGED  TO  0A9S. 


r • • i r»  * * • 


O ■>  j 1 o '-6  o ? o rb  I O o 


CCM=>>JT5  EIGENVALUE  - 5 N 0 c I G IN  V EC  TORS  OR  3 


< = 1 

L=M* (U-l > 

do- too  i-i,j 

L = L * 1 
< = <*1 

IOC  3 ( K)  =3  (<_  I 


THE  lAT^IX  9 IS  a REAL  symmetric  matrix. 
MV  = 0 

0Att~r TSE N (0,  X,M > MV )■ 


ARE  RR2HULTI=LIED  BY  THE  ASSOCIATED  EIGENVECTORS. 


L = 0 

00"  t~10~  J^ItM- 

L = L + J 

~rt o- - vti-j  r =i . 'TV — s m 

K=0 

BC-ttE  i=l ,M 

00  115  1=1, M 

K-K-^l 

HE  3 (<)  =X  (<)  * XL  ( J ) 


fq=>M  (9»*(-l/2))  PRIME  * A * <3**(-l/2)) 


DO  120  1=1, M 
DO  120  J = 1 , M 

M !->»»(  1-1  y 

1=“* (J-l ) +1 

X-(  h 0 = 0 » 0 

30  12  0 < = 1 , M 

N1=N1+-1 

N 2 = N 2 + 1 

— 1-2 »T- X- f L V-X  YL  h+  f M2  ) 

L = 0 

- 00-130  0 = 1,  ** 

DO  130  1=1, J 

M2=M* ( J- 1 ) 

tr=trH 

A (L) =0  . 0 

9013:3  < = 1,M  — 

N 1 = N 1 + M 
H 2=M  2 + 1 

130  A (L)  =A  (L)  +X('I1  1*9  <N2> 

C OO^R'JT'  EIGENVALUES  AND  EIGENVECTORS  Cc  A 

CA._l  EIGEN  (A,X,M,MV) 


140  Xl<I)=A(L) 


C OOMPUTr  TH"  NOPNALIZED  EIGENVECTORS 

C 

no  iso  : = i,m 

N2=0 

00  ISO  J=i,M 
N1=I-M 
L = M* (j-l ) +1 
A (L) =0.0 
00  150  K = 1 , M 
N 1 =N 1 + M 
N2=N2+1 

150  A(U=A<L>+0(N1)*X(N2) 

L = C 
K = 0 

00  100  J=1,M 
SUMV=0.fl 
no  170  1=1,4 

L =L+  1 

170  SUMV=SU-1V  + a (L  ) *A  (L) 

17  s SUMV=  SORT ( SU  MV) 

oo  lao  i=i, m 

K =K+  1 

180  X IK)  =A  (< I /SUMV 
RETURN 
c NO 

I 


o c 6 o Im.j  6 o t I I I I o (bu  I 1 I I ! c»o  t)  o o 


r 


THE  DOUBLE  DC?ECISIOK  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 
CONTAIN  00 U BL  E PRECISION  FO°TPA  N FUNCTIONS.  SQPT  IN  STATEMENTS 
43,  71,  4NC-7S“MUST  BE' CHANGED  TO  OSQRT . ABS  IN  STATEMENT 

52  MUST  BE  CHANGED  TO  DABS. 


GENERATE  IDENTITY  MATRIX 


IF(MV-l)  10,25,10 

— ID  ID--N 

DO  20  J=1,N 

10— ID^N — 

DO  20  1=1, N 

IJ=IO*I 

R ( IJ  ) = 0.  0 

lc  > g-O-,15 , 20 

15  R ( I J ) = 1.  0 

— -eOHT  I HUE — 

eOHPU-T E-  INITTAt-ftNO-FI NAt^NORHS  (ANORM  ANQ— ANORHXD 


— g-"  A NORM=  0 . 0 

DO  35  1=1, N 

- - DO -35 — U=I,H 

IF(I-J)  30,35,30 

30  - 1 A=I  

ANOPM=A.ICRM  + A (IA)*A(IA) 

—3=  COMT  IN  US 

IF { A NO  PH ) 165,165,40 

—4  0— *NO‘RM=lrr-4-tV^SO‘RT-(-ANORH) 

AN^MX=  ANORM*! . 0E-6/FUO AT (N) 


INITIALIZE  INO ICA  TOPS  ANO  COMMUTE  THRESHOLD,  THR 


I n o = n 

THR=ANORM 

4"  THPsTHO/FLCAT  (N) 

50- t-t 

^ Mr  L+ 1 


7 

I 


> 


COMPUTE  SIN  ANC  COS 


60  MQ=(M*M-M>/2 

LO=<L*L-LT-A2 

LM=L«-MQ 

6~2 — I~F  ( - ABS  ( A — -T-HR-) — IS  OtBEt&E 

65  INO=l 

LL-tM-tO 

mm=m+mo 

X = 0.  5*  H ( t L >—  A fMM-H 

6 A Y = - A { L M) / SORT  (A <LM) *A (LM)+X*X) 

t-F-FX->— TO-,  ?5t-?-5- 

73  Yr-y 

75-  SINX  = Y/-SC;JT(2.0*(1.  0+  ( - SORT  < 1 . 0-Y*  Y T T-H — - 

sinx2=s:nx*sinx 

7^-OOS-X— SORT-t  1.  3-5-1  NX  2 > 


BE 


* « < H J 


m copy 
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C05X2=C0SX*C0SX 
SINCS  - 5INX*COSX 


r 

r 


I LQ=  N* ( L - 1 » 

IMQ  = N*  (M-l  ) 

00  125  1 = 1, N 
10= ( I* I- I » /2 
IF(I-L)  80,115,80 
80  TF(I-M)  85,115,90 
8 r’  IM  = I«-MO 
GO  TO  95 
90  T^rM+TQ 

95  IF(I-L)  100,105,105 
100  T L = I + L Q 
GO  TO  110 
105  I L=L  +1 Q 

110  X=A(IL >*COSX-A  (IH)*SIMX 

A(TM)=A(IL)*STNX+A(IM)*COSX 
A (IH=X 

llc  IF(HV-l)  120,125,120 
120  IL<R=ILQ«-I 
IMR=  IMQ-*- 1 

X=R(ILP)*COSX-R(IHR)*SINX 
R(IMR)=R(ILR>*SINX4-R(IMR)*C0SX 
R ( TLR ) =X 
1 2 r CONTINUE 

X = 2.0*A(I_M)*SINCS 
Y=A(LL)*COSX2+A(MM)*SINX2-X 
X= A ( LL )*SINX2+A(MM)*C0SX2+X 

A(LM)=(A(LL)-A(HM) >*SINCS+A(LM)* (C0SX2-SINX2) 
A (LL)  = Y 
A ( MM ) = X 
C 

C TESTS  FOR  COMPLETION 

C 

C TEST  FOR  M = LAST  COLUMN 

C 

130  TF(M-N)  135,1*0,135 
135  M=M+i 

GO  TO  GO 
C 

C TEST  FOR  L = SECOND  FROM  LAST  COLUMN 

C 

140  IF(L-(N-1)>  1*5,150,145 
145  L=L4-1 
GO  TO 

130  TF(INO-l)  160,155,160 
155  INO=  0 

GO  TC  5] 

C 

C COMPARE  THRESHOLD  WITH  FINAL  NORM 

C 

160  IF (THR-ANRMX)  165,166,46 


C SORT  EIGENVALUES  ANO  EIGENVECTORS 


c 


16"  IQ=-N 

no  185  1=1,  n 

T 0=  IQ-f  N 

LL  = I-*-(  I*I-I)/2 
JQ=N*{ 1-2) 

00  185  J=  I ,N 
JQ=JQ*N 

MM-J+( J‘J-J)/2 
rf(A(LL)-A  (MM ) ) 170, 18?, 185 
170  X=A ( LL ) 

A ( LL ) = A ( H M ) 

A ( MM ) = X 

IP(MV-l)  l-7?, 185,175 

1 7 c no  180  <= 1 ,N 

ILP= IQ  +< 

TMP= JQ ♦< 

X = P(IL=>) 

R(IL9) =R  < I HR > 

180  P ( IMP ) =X 
18?  CONTIMU" 

»rTUPN 

1ND 


APPENDIX  H 


SUBROUTINE  CHECK 


SUBROUTINE  CH'CK , RETURNS  (R1,R2> 

REAL  LX,LXO,LY,LYB,LXY,NX,NXB,NY,NYB 

COMMON/CNTL/INTER,PREV, VAL,  1ST ART, K, M A X , NRUN,  ICONV  , T IME  1 , T IME2 
CONMON/SOLN/AA  (2  4 01)  ,6B(2401> ,EIG<49> , EIGV<2401) ,MOOE 
CO^MON/LCAD/L X ,LXP,LY ,LYB,LXY,NX,NXO,NY,NYQ 
LOGICAL  INTER 
TOLR=. 002? 

K2=K*< 

IF(  ASS  CIG  (1)  > .LT  . 1.  E-40  >EIG<1)  = 1. E-4  0 

~ ig ( i ) =i./-:ig  id 

V AL=EIG(1) 

HOPE  =1 

00  1 1=2, K2 

IF (AOS ( :IG (I> ) .LT. l.E-40 )EIG  <I>  = 1. E-40 
E IG  ( I) = 1 . / ~ IG  (I) 

IF  (ABS ( IIG (I) > .Gr. A8S< VAL) > GO  TO  1 
V AL  = EI G ( I > 

NODE  = 1 

1 CONTINU' 

IF (1ST  ART . GT . 0 >G0  TO  10 

IF (AOS ( V A L ) .EQ.ABS(PREV) ) I C 0 N V = 1 

0 IFF  = A OS  ( ( ABS  ( VA  L ) -A  BS  (PREV)  > /VAL) 

IF (OIFF. LT .TOLR> IC0NV=1 
IF(.NOT. INTER)CALL  PRINT 
°R?V=VAL 
< = K+  1 

IF  (IC0NV.EQ«1)  RETURN  R2 
IF (K.LE. MAXJRETURN  R1 
WRITE(G,900)OIFF 
RETURN  <2 
10  ISTART  = 0 

IF (. NOT. INTER) CALL  PRINT 
pp  = V = V AL 
K = K+1 
RETURN  ’1 

900  FORMAT  (/10X,-,2H**  WARNING*  RELATIVE  CONVERGENCE  TOLERANCE, 

1 /13X,22HOF  .002?  NOT  SATISFIED, 

2 / l 0 X , 1 * TOLEcANCE  OF  ,F7.4,9H  ACHIEVED/) 
cNr 


) 


appendix  I 
SUBROUTINE  OUTPUT 


* 


> 


i « 
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J 


SUBROUTINE  OUTPUT  , RE  TURNS (R1 ,R2) 

’24L  LX,LX3,LY,LY0,LXY,NX,NX9,NY,NY3 

COMMON /M  A Ttr/S  C2T , T (2  ) , XNU  ( 21  , XLMO*  (2\ , BE T A f2  ) , S IGN  (2> 
C03M0N/CCRS/GX,GY,TC,A,B 

CtrwrO'r/^Nrc/ruYr^TF><?i-V-rV'ACvr2Trrfrr,-)crM'^'«%-‘TPt;?rrIcc'TU-rTrn£-ir-,-T'T'*r2- 

C0MM0N/30LN/A A (2401)  , BB  < 2401  > ,EIG(49),EIGV<2401>,MOOE 
COMMON/LCAO/L  X ,LX?,L Y,LY9,LXY,NX,NX3,NY ,NY9 
LOGICAL  INTER 

'TATA  Y?.G/1PT/ 

NRIJN  = N PUN- 1 

K-r-K’-l 

ISTART  =3 

W?-;ET~5-t-6  , »9fr>X 

<2=<*< 

— X=VAL*NX^ — 

Y=VAL*NY 

XY  = V'AL*LXY 

XB=V AL’NXB 

Y 9-V  At*NY8 

WRITE( o, 91 3) 

X°TT?(  6 , BTrEhX-,  XY~,  Y 

WRITE(6,91?)X3,Y3 

f F t , NO  T .1  NTER-) — DO — TU — 1~0~0 

WR ITE ( 6 , 995) 

READ (5, 3 97)  A NS 

IF(ANS.'1E.YE5)RETURN  R 2 

13  0 — CONTINUE' — 

WRITE ( 6,9  0 0) 

W-RUTE  16 , 98"8T  f = IS“YUUT  » J"J=  t,K-2T 

J -K2  * ( MOE-1)  +1 

OJ=i<2*MOOE 

WRITE ( 6 » 92  0 ) 

WRITEf6,919) 

IF(K.LT.10)WRITE(6,940)K,K,K 

W RETtM  6t9  25  YfE'IG  AYrt  X-rT^UvUJ  > 

IF(NRUN.LE.O) RETURN  R2 

--RETURN-  Eir 

?90  FORMAT  (13X  ,23H*****  SOLUTION  FOR  M=N=,I2,5H  •***) 

39?  FO°HAT (//10X, 40HPOINT  EIGENVALUES  AND -MODE- SHAPES  tY/N)?»— 

3R7  FOcM AT  ( A5 ) 

930 F-CPXAT-E/UO  Xy-29H*~»-*- **-*■*-*  *-  EI-GE-NVfrtUES— *-»•♦  ♦ »**  ♦ > 

905  c0=mat  (3X,^?12.5) 

910  FORMAT  (/ 1 0 X , 30  H***  ***•-  CRITICAL  LOADS  ***•**-*) 

915  FORMAT  (1GX  ,4HNX  = , FI  6 . 4 , / 10  X , 4HNXY  = , F 1 S . 4 , / 10  X , 4HN  Y =,F16.4) 

917  cOPMAT  (13X  ,4HNXB- ,Fie.  4 , HH  *<T-2*-Y/3>, 

1 /13X  ,4HNYB=  , F15.  4 , HH  *(1-2*X/A)> 

— 9t9 F0oHjT-ft7x-,i’-HiT-eRI-TICAt-LCA0S-) 

920  RQFM AT (/ 10 X , 30** ****** **  MOOE  SHAPE  **«*♦***») 

9 25  FORMAT  f=  Y,  4F±B. 10  ) 

9-»  0 FORMAT (/ 10 X,1 5H( Wll, W12, ...  ,W1, 1 i, 14H,W21 ,W22, ...»  W, 211, 1H> /) 

ENTRY  ECHO  

WRITE!  6,  94  3) 

WRIT E t 6,-949) 

WRITE  <6, 95  0>£  (1)  , XNU  (U  , T (1 ) 

WRIT?  ( 6. 95  t)E  <2>  , XNU  <2  ) , T ( 2) 

W R IT  E ( 6,  947) 

W»IT-(6,952>GX,GY,TC 
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WRITE  < 6,  9<*6) 

WRIT"  (6*  95  3)A 
WRITE  < 6, 954)  B 
WRITE  ( 6,  945) 

WRITF(6,955)NX,LXY,NY 

WRITE(6,957)NXB,NYB 

IF  C. NOT. INTER) WRITE ( 6, 960) 

945  FOPMAT </5X ,12HLOAO  FACTORS) 

9-*  6 FORMAT (/5X ,19HPLANF0RM  DIMENSIONS) 

9-*  7 FORMAT  (/5X  ,13HC0RE  PARAMETERS) 

9^8  FORMAT  ( 1 H 1 » /5 X t2 1 HFACE  SHEET  PARAMETERS) 

949  FORMAT </23X,lHE,10X,2HNU,8X,iHT) 

950  FORMAT  (3X,7HFACE  1 1 , 2X , E14 . 4 , 2F 1 0 . 5 ) 

951  FORMAT <5X,7HFAC£  2 : * 2X , E 14 . 4 , 2F1 0 . 5 ) 

952  FORMAT (/22X,2HGX , 12X,2HGY, 10 X,2HTC , 

1 / 5X,7HCOR£  t ,2X,2E14.4»F10.5) 

951  FORM  AT  (/5X  ,7H  A =,F10.4> 

954  FORMAT ( 5X,7H  B =,F10.4) 

955  FORMAT {/5X,3HNX=,FE. 2 , 3 X , 4HN X Y= , F5 . 2 , 3 X , 3HNY=,F5.2/) 

95  7 FOFMAT  (3X,4HNXB=,F5.  2 , 14 X , 4HN Y R= , F 5 . 2/ ) 

960  FORM  AT  (/5X ,5HTERMS,7X, 2 HN X, 7 X , 3HNX B , 8 X , 3HN X Y , 7 X , 

1 2HNY ,8  X, 3HNYB , 3X , 2 OH ASSE M < SEC ) E IGE N < SE C ) ) 

RETURN 
ENTRY  PRINT 
X = V A L*  NX 
Y=VAL*NY 
X Y = V AL  *L  X Y 
X B = V AL  *N  X B 
YP=V AL’NYB 

WRITE(6,970)K,X,XB,XY,Y,YB,TIME1,TIME2 

970  FORMAT  (6X,I2i2X*TF10.3l 
RETURN 
END 
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APPENDIX  J 


SUBROUTINE  BENDI 


SUBROUTINE  BEND1 

REAL  LX,LXB,LY,LYB,LXY,NX,NXB,NY,NY3 
C0MMCN/C0RE/GX,GY,TC,A,9 

COMMON/: NT L/INTER,PREV,VAL, I START, K,MAX,NRUN, ICONV , TIME  1 , TIME2 

COMMON/SOtN/4  A (2401)  , (30(24  01)  , EIG(49),EIGV(2401)  ,MOOE 

C0MM0N/LCA0/LX,LXB,LY,LYB,LXY,NX,NX9,NY,NYB 

LOGICAL  MREO, NSEQ,MR00 D.NSOOO,  INTER 

OAT  ft  PI/3. 1*13926535898/ 

00  100  “1  = 1, K 
00  100  N=  1 ,K 
no  100  IR=1,K 
M RE  Q = ( M. EQ. IR ) 

I CHK  = ( M+ I R ) /2 
XCMK=. B’FLOAT (M+IP) 

MROOO=  (A  OS (FLOAT  (ICHK)  -XCHK)  ,GT.  .10) 

no  ioo  i s = i , k 

NSCQ=  ( N. EO.IS) 

ICHK= ( N+ IS ) /2 
XCHK  = , 5*FL  OAT  ( N+  I S ) 

NSOOO=  (A BS  (FLOAT  (ICHK) -XCHK) ,GT. .10) 

ICHK= ( IR-1)*K-HS-1 
JJ=ICHK»K*K+ (M-l ) *K+N 

IF ( ( .NOT ,MREQ) .AND. ( .NOT ,NSEQ) )GO  TO  100 
IF(MREQ. ANO.NSEQ) BR( J J ) = BB  ( J J ) - . 12  5 *P  I *P  I * 

1 (M*M*B*LX8/A*N*N*A*LYB/B>*4./ (A*8) 

IF  (MREG.  ANO.NSOnO)  BB  ( JJ)  =B8(  JJH-2.  *M*M*N* I S*LXB 
1 *B/(AMN*N-IS*IS)**2)*4./(A*B> 

IF (M ROOT. A NO.  N SE Q ) BB ( J J ) =B8 ( J J ) ♦ 2. *N*N* M* IR*L YB 
1 *A/(BMM*M-IR*IR)**2)*4./(A*B> 

100  CONTINUI 
RETURN 
END 
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APPENDIX  K 


SUBROUTINE  ALTER 
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‘.1  1 'Hi.1  -V 


SUBROUTINE  ALT6R, ecTURNS (R1  ,R2,R3) 

RCA L LX,LXn,LY,LYe,LXY,NX,NXB,NY,NY9 

COMMON/MATL/E (2)  , T (2 ) , XNU(2> ,XLMOA (2) , BET  A (2) .SIGN (2) 
COMMON/OORS/GX  ,GY  , TC  , A ,B 

COMMON /O NT L/INTSR,PRcV,  VAL,  I S T ART , K , M A x , N»UN , I C ON V , T I ME  1 .TIME 2 
CCMMON/3CLN/AA(2401),Bn(24Gl>,FIG<49),rIG7(2401> ,MOOF 
COMMON/LCAO/L  X, l XB.LY.LYO.LXY, NX, NX9,  NY, NY  IL- 
LOGICAL INTER 
0 T MENS  I IN  VAR(17) 

OATA  YtS/lHY/ 


OATA 

7 AR/4HE 1 

, 4 HL2 

, 4HNU1 

, CHNU? 

1 

4 HT  1 

, 4HT2 

, 

2 

i*  HGX 

, 4 HGY 

, 4 H T C 

3 

4 HA 

, 4 H 3 

» 

u 

AHNX 

, 4HNX  Y 

, 4HN  Y 

* 

5 

4 HNX  9 

, 4HNYB 

, 4HENO 

/ 

K0UNT=0 

WPITE(G,900) 

WRITE <6, 910) 

PEAO  (5  , 120  ) ANS 

IF  (ANS. EQ.YSS)  GO  TO  10 

WRTTE(6,960) 

RFAO(5,920>  ANS 
IF (ANS.  \Q. YES)  RETURN  PI 
PE  TURN  R3 
10  CONTINUE 

WRIT" (6,930) 

15  CONTINUE 

WRITE(6,940) 

RPAnf-  ,030)7 

IF (V . EO  V A R (1 T ) ) GO  TO  800 
OO  20  1=1,17 

IF (V.NE. VAR (I ) ) GO  TO  20 
J = T 

GO  TO  2> 

20  CONTINUE 

W R I T E ( 6 , 9 7 0 ) 7 
K OUN  T=  KOUN  T-*- 1 
IF  (KOUNT.GT. 3)  RETURN  R3 
WRITE (6,975) 


GO  TO  1) 


23 

WRTTE 

(6,990)  7 

° E AO  ( 

5 ,* ) 7AL 

GO  TO 

(101,102 

1 

109,110 

101 

E (1)  = 

7 AL 

GO  TO 

1 ' 

10? 

E (2)  = 

7 AL 

V, 

GO  TO 

1 •> 

10  3 

XNU( 1) =/AL 

GO  TO 

1 - 

4 

10- 

X NU  ( 2 

) =7AL 

.* 

GO  TO 

1 ’> 

, 

105 

T (1)  = 

7 AL 

v> 

GO  TO 

1 • 

105 

T (2)  = 

7 AL 

GO  TO 

i> 

65 


- -■  :■  - 


ZA 


GX=VAL 
GO  TO  l' 

GY=VAL 
GO  TO  15 
T C = V AL 
GO  TO  15 
A = VAL 
GO  TO  15 

0 = VAL 

GO  TO  15 
NX-  VAL 
GO  TO  15 
LXY=  VAL 
GO  TO  15 
NY=  VAL 
GO  TO  15 
NXR=  VAL 
GO  TO  15 
NYR=  VAL 
GO  TO  15 

continue 

SET  RESTART  PARAMETERS 
NRUN=0 

1 CON V = 0 
P»EV=0. 

CALL  RESET 
WRITE (6,983) 

WRITE(6,900) 

RETURN  R2 

FORMATS 

FORMAT </10X,  11HBEGIN  A L T ER , / 1 0 X , 17 HMOO I FY  DATACY/N)?) 
R0  CM  A T <AL) 

FORMAT (10X,20HGIVE  ^PARAMETER  NAME) 

FORMAT  ( 15 X » 1 H* ) 

FORMAT ( A 5 ) 

FORMAT ( / 1 0 X , 2 3 HB EG IN  NEW  PROBLEM (Y/N ) ? ) 

FORMAT (/10X,14HINVALID  OATA  »,A5/) 

RQRMAT  ( 1 0 X » 1 7 HRE  T RY  OP  TYPE  END  /IX) 

FOPMATf  10X,15HALTER  COMPLETED) 

ROPMAT (16X ,A5 , 3H=  ) 

END 


APPENDIX  L 
SUBROUTINE  TITLE 


SUBROUTINE  TITLE  (I) 

1 r 0 c M A T ( / / ) 

2 R0RMAT(19X,43H  ****************************  *♦**♦****♦♦***♦) 

3 POpMAT  (18X ,A3H*  *> 

£.  FORMAT (18X  ,L3H*  GENERAL  STABILITY  ANALYSIS  OF  ♦> 

5 FO?M AT  (l 9X ,43H*  FLAT,  RECTANGULAR  SANOWICH  PANELS  *) 

6 FORMAT C 1 a X ,A3H*  REVISEO  10/15/76  ♦) 

7 FORMAT (18X,43H*  UNIVERSITY  OF  DAYTON  RESEARCH  INSTITUTE  *> 

8 FORMAT  (13X  ,43H*  DAYTON,  OHIO  *) 


9 FOPMAT(lHl) 

IF  (I.LE.O)  WRITE (6,9) 
WRITE (6, 1) 

IF  (I.GT.l)  RETURN 
WRITE(6,2) 

WRITE(6, 3) 

WRITE (6, A ) 

WRITE (6,5) 

WRITE(6,5) 

WRITE( 6, 3) 

WRITF ( 6, 7 ) 

WRITE(6,8  > 

WRITE (6, 3) 

WRITE(6, 2) 

WRIT; (6, 1 ) 

RETURN 

END 


